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Motivated by recent progress in epitaxial growth of proximity structures of s-wave superconductors 
(S) and spin-active materials (M), we show that the periodic structure of S and M can behave 
effectively as a superconductor with pairs of point nodes, near which the low energy excitations are 
Weyl fermions. A simple toy model, where M is described by a Kronig-Penney potential with both 
spin-orbit coupling and exchange field, is proposed and solved to obtain the phase diagram of the 
nodal structure, the spin texture of the Weyl fermions, as well as the zero energy surface states in 
the form of open Fermi lines (“Fermi arcs”). Going beyond the simple model, a lattice model with 
alternating layers of S and magnetic Z 2 topological insulators (M) is solved. The calculated spectrum 
confirms previous prediction of Weyl nodes based on tunneling Hamiltonian of Dirac electrons. Our 
results provide further evidence that periodic structures of S and M are well suited for engineering 
gapless topological superconductors. 


The time-honored recipe for discovering new supercon¬ 
ductors with interesting pairing symmetries or topolog¬ 
ical properties is via the synthesis of new materials. In 
recent years, an alternative approach has been advocated 
and gained experimental success. It is based on making 
proximity structures of s-wave superconductor (referred 
to as S hereafter) and spin active materials (M) such as 
semiconductors, topological insulators^’^ or ferromagnets 
with spin-orbit and/or exchange coupling. With proper 
design, the proximity structure can behave effectively as 
a superconductor with the desired symmetry or topol¬ 
ogy, at energies below the bulk superconducting gap of 
S. For example, Fu and Kane^ showed that the interface 
of a three-dimensional topological band insulator (TI) 
and an s-wave superconductor is analogous to a spin¬ 
less Px + ipy superconductor that hosts Major ana zero 
modes at vortex cores. Similar states also arise in prox¬ 
imity structures of S and two dimensional electron gas 
with Rashba spin-orbit coupling and Zeeman splitting 
either due to a nearby ferromagnetic insulator or an ex¬ 
ternal magnetic field^“^. In one dimension, e.g., a semi¬ 
conductor nanowire^ or a chain of ferromagnetic atoms 
deposited on a superconductor^^, Majorana zero modes 
form at the sample edges. While it remains a challenge to 
fabricate and control the interface properties, or detect 
the unequivocal experimental signatures of these states, 
significant experimental progress has been made in recent 
years (for an overview, see Refs. 11 and 12). 

In this paper, we explore the possibility of realizing 
gapless topological superconductivity in S-M proximity 
structures. Specifically, we focus on superconducting 
states with topologically protected point nodes, i.e., the 
analogs of the A phase of superfluid ^He^^ and Weyl 
semimetals^^”^^. Such a state was predicted to appear in 
the superlattice structures of superconductors and mag¬ 
netic topological insulators by Meng and Balents, and re¬ 
ferred to as “Weyl superconductors”^^. The elegant anal¬ 
ysis of Ref. 17 is based on an effective Hamiltonian de¬ 
scribing the tunneling of the helical Dirac electrons (the 
surface states of TI) across the layers of TI and S, where 
the presence of the superconducting pairing potential and 


Zeeman field provide a mass to the Dirac electrons. The 
proposal of Meng and Balents can be viewed as the gen¬ 
eralization of the earlier work on Weyl semimetal in the 
multilayer structures of trivial insulators and magnetic 
topological insulators^^. One is then led to the following 
questions: is it feasible to realize Weyl superconductors 
using materials other than topological insulators; is heli¬ 
cal Dirac electron essential? 

We answer these questions by considering a simple, 
idealized model of S-M superlattice. Here M stands for a 
general spin active material with both spin-orbit coupling 
and exchange splitting (i.e., both the time reversal and 
spatial inversion symmetry are broken). We assume, as 
in Ref. 17, that the M layer is sufficiently thin so that 
the suppression of superconductivity is not significant 
and the superconducting phase coherence is maintained 
across the M layers. This motivates us to approximate 
the M layers as delta function spin-active potentials, sim¬ 
ilar to the well known Kronig-Penney model. The band 
structure of this model is solved to illustrate the recon¬ 
struction of the low energy spectrum due to the peri¬ 
odic spin-active potential. We identify parameter regimes 
where the spectrum has one pair or two pairs of Weyl 
nodes. We discuss the low energy effective Hamiltonians 
near the Weyl nodes, and the Fermi line (usually referred 
to as “Fermi arc”^^) zero energy surface states that mani¬ 
fest the nontrivial topological properties of the supercon¬ 
ducting state. This model clearly illustrates that neither 
TI nor helical Dirac electron is necessary for Weyl nodes 
to appear in S-M superlattices. 

This simple model can be straightforwardly general¬ 
ized to treat finite thickness of the M layer. It however 
neglects many microscopic details of the specific materi¬ 
als and the S-TI interface. For the purpose of providing 
design parameters for Weyl superconductors, it is also 
desirable to describe the periodic structure by tight bind¬ 
ing models defined on discrete lattices. We present such 
a lattice model for the S-TI superlattice, in which each 
unit cell consists of a few layers of S and another few 
layers of magnetized TI with a tunable hopping matrix 
describing the coupling between the two materials. We 


2 


outline a procedure to compute the energy spectrum of 
the superlattice, and assess the requirements to realize 
Weyl superconductors with one pair of nodes and two 
pairs of nodes respectively. Along the way, we briefly 
review the properties of a single S-TI interface^, and dis¬ 
cuss the relation between the Andreev bound states at a 
single interface^^ and the spectrum of multilayer systems. 


I. A SIMPLE MODEL FOR S-M 
SUPERLATTICE 

We consider a periodic layered structure of an s-wave 
superconductor (S) and a spin active material (M) ex¬ 
tending in the 2 : direction as schematically shown in 
Fig. 1. As far as the low energy excitations are concerned 
(relative to the bulk superconducting gap A), M amounts 
to a periodic spin-active potential V(z^d) = V(z), where 
the hat denotes matrices in spin space. To preserve the 
superconductivity throughout the whole structure, the 
M layers should not be too thick so we assume the M 
layer thickness is much smaller than the period d. In this 
limit, V can be modeled by Kronig-Penney potential of 
the form 

r(k||, = d ^ (5{^ - nd) [Fo + ^ (k|| )ay (1) 

n i 

Here d’s are the Pauli matrices in the spin space, ky = 
{kx^ky) is the transverse momentum which is conserved 
due to translational invariance on the xy plane. Note 
that the material details of M do not enter in this de¬ 
scription, they are encoded in Vq and Vi which are cho¬ 
sen to reproduce the scattering matrix of electrons by 
M. The superlattice is then described by the following 
Hamiltonian in the particle-hole space, 

ho{ki\,z) ^ Ai(T2 1 ^2) 

-Az( 72 -/io(-k||,2;) J ’ 

with /io(k||,2:) = (k| — d^)l2me — ^ + V{k\\,z), and 
the check denotes a matrix in the particle-hole and spin 
space. Note that we have assumed A to be homogeneous, 
and for ^ = nd, i.e. insides M, the potential V domi¬ 
nates all other terms in the Hamiltonian. For simplicity, 
we shall put Vb = 0- To model M with both spin-orbit 
coupling and exchange splitting, we assume that V takes 
the form^^ 


^(k|| 


Vi{kii)ai = Vsoi-kyOx + k^Oy) + (3) 

i 

where Vso is the strength of Rashba spin-orbit coupling, 
and Vz is the Zeeman (exchange) field along the z direc¬ 
tion. Using the periodicity of the Hamiltonian, 1-L{z) = 
ld{z -h d), the band structure can be easily obtained by 
the expansion of the wave function via Bloch’s theorem 

y, z) = Y. (4) 

G 


where G is reciprocal lattice vector G = 7 x 1211 jd and m G 
Z. The generalized Bogliubov-de Gennes (BdG) equation 
then becomes 


^ [%(k,G)(5G,G' + V(k||)]l>k.G' = (5) 

G' 


Here we have separated the “unperturbed” Hamiltonian 


no{k,G) 


^{k\\,kz + G) Az (72 

—Aia2 —^(k||, + G) _ ’ 


(6) 


with ^(k||, + G) = (k| + + G)‘^)/2me — y, and the 

spin-active “perturbation” 

V(k||) = Vsoi-kya^ + (g) k^ay] + O cr^ (7) 


with Tz being the Pauli matrix in the particle-hole space 
(we will drop the hat for a when there is no ambiguity). 



FIG. 1. Schematic of the S-M periodic structure. 


The infinite dimensional matrix equation in Eq. 5 can 
be solved numerically by a truncation, keeping only \G\ 
up to some large enough value of Nir /d, followed by diag- 
onalization to yield the band dispersion E\^. This trunca¬ 
tion is physically equivalent to introducing a small width 
to the M layers. Of course one has to check that the 
low energy spectrum does not depend on N. This model 
can easily be generalized to the case of M layers with fi¬ 
nite thickness. For example, superlattice unit cell can be 
modeled in a way that the region z G [0,di] is occupied 
with S (where A is constant, V is zero) and 2 ; G [di, d] is 
occupied with M (where A vanishes but V is constant). 
In this case, both A and V have off-diagonal matrix el¬ 
ements in G space and the resulting BdG equation is 
slightly more complicated than Eq. 5. 

We are particularly interested in the zero energy solu¬ 
tions of Eq. 5. Eor this purpose, it is useful to introduce 
0k = which can be shown to satisfy the follow¬ 

ing equation; 

ik.s0k= [l-^(S-?7o(k,G))-iV(k||)]0k = O. (8) 

G 

The matrix inverse is the bare Green function of the bulk 
superconductor and can be computed analytically. Ex¬ 
istence of the zero energy solutions {E = 0) at isolated 
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nodal points on the kz axis with ky = 0 is equivalent to 
the existence of a non-trivial solution of Eq. 8 which can 
be expressed as det^dk^o = 0- This equation can be fur¬ 
ther manipulated analytically to give the following simple 
equation as the condition of the zero energy solutions; 

g^ + f-I = Q, (9) 

where g = V^Y,m ^rn/ (C™ +A2), f = v^Y.m A/(Cm +A2) 
and ^rn = {kz -^27rmld)‘^ j2me —/i. This offers a fast way 
to scan for Weyl nodes in the parameter space since no 
numerical diagonalization is required. 

Before we present any numerical results, it is worth¬ 
while to develop a qualitative picture of the low energy 
excitations in such S-M superlattices. Consider the nor¬ 
mal state dispersion (turning off superconductivity by 
setting A to zero) in the absence of M. For ky = 0, 
the low energy excitations are located at large momenta 
around kz ^ :^kF. In the presence of a weak periodic 
potential E, the spectrum of the superlattice structure 
can be obtained by folding the free dispersion into the 
first Brillouin zone kz G [—Tr/d, Tr/d] which gives a set 
of Bloch bands, {kz + G)^/2me, all being restricted to 
small momenta since n/d kp- The spectrum acquires 
a gap as A is turned on. A finite Zeeman field will split 
the Andreev bound states formed below the supercon¬ 
ducting gap A, and push one of the branch towards the 
zero energy. A non-zero spin-orbit coupling can in prin¬ 
ciple endow a topologically nontrivial spin structure to 
these zero energy states. The S-M proximity structure 
considered here differs from the well studied system of 
semiconductor nanowires^ in its dimensionality. In three 
dimensions, linear dispersion in the vicinity of the nodes 
is known to be described by the Weyl Hamiltonian^^dQ 


II. WEYL NODES 

Two representative examples of the low energy spec¬ 
trum of the S-M superlattice are shown in Fig. 2. We 
will explicitly show that they correspond to one pair and 
two pairs of Weyl nodes, respectively. In both cases, we 
observe a linear energy-momentum dispersion near iso¬ 
lated points {k^} located on the kz axis. At any of these 
diabolical points, the zero energy state is doubly degen¬ 
erate. Let us label these two degenerate states as |T+) 
and |T_), or |±) for short. The low energy physics near 
the node is described by an effective Hamiltonian which 
is a 2 X 2 matrix in the Hilbert space spanned by |±) with 
the general form^^ 

1 ^ 

^(q) = = F q = k - k°. (10) 

hi=i 

Here fi are the pseudospin Pauli matrix in the ± space. 
The low energy Bogoliubov quasiparticles thus resemble 
massless chiral fermions (Weyl fermions). The chirality 


refers to the locking of the pseudospin with respect to the 
momentum direction as described by h(q). The direction 
of h, /i(q), constitutes a mapping from a sphere in the 
momentum space enclosing the Weyl nodes to an unit 
sphere in pseudospin space. The topological invariant 
for this mapping is the winding number^^ 

^ = T Cyfc f d^kh ■ {dji X djh) (11) 

where the integration is over a closed volume containing 
q = 0^ and dj = d/dqj. In the simplest example h = ±q, 
= ±1; and for Vij = Xidij^ = n^sign(Ai). Since 
resembles to a spin-1/2 particle in magnetic field, one 
can define Berry connection and the corresponding fiux 
densityThen a Weyl node corresponds to a magnetic 
monopole with charge in the momentum space. For 
periodic systems, the net flux through the Brillouin zone 
(BZ) and the net magnetic charge inside the BZ must be 
zero^^. Therefore Weyl nodes always appear in pairs of 
opposite charge They are well separated in k space 
and topologically stable^^. 

E/A 

;a.oo 


E/A 



TV 


FIG. 2. Weyl nodes in S-M superlattices. Only the low¬ 
est eigen energy E{kx,ky = 0,kz) is shown. Blue line is its 
dispersion along the kz axis inside the first Brillouin zone 
kz G [—7r/d, 7r/d]. Left panel: single pair of Weyl nodes. The 
range of kx is kxjkp G [—.05, .05]. d = IOtt/Zcf, Vz = 2.5A. 
Right panel: two pairs of Weyl node for kxjkp G [—.015, .015], 
d = 9.57r/A:F, Vz = A. Here d is the super lattice period, Vz is 
the strength of Zeeman field and Vso is the spin-orbit coupling. 
We choose A = 0.05/i and Vso = 5A. 

It is important to conduct the search for zero energy 
states in the entire momentum space. Fig. 3 shows 
the energy-momentum dispersion along the kx axis for 
fixed kz = k^ for various values of Vgo- For very small 
spin-orbit coupling, the slope around kx = 0 is vanish¬ 
ingly small, and there are many other low lying states 
at larger values of kx that are sufficiently close to zero 
energy. These are Andreev bound states (ABS) ubiq¬ 
uitously found in superconducting proximity structures. 
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For example, for Vgo = 0, the structure reduces to 
a superconductor-ferromagnet (S-F) superlattice. It is 
then expected from semiclassical consideration that a se¬ 
ries of ABS with finite kx will be formed between two 
adjacent F layers. The slope near kx = 0 increases with 
Vso- At the same time, other ABS at finite kx are increas¬ 
ingly gapped out. For large enough Vgo, the Weyl nodes 
on the kz axis are the only zero energy states. In other 
words, spin-orbit coupling is crucial for the S-M super¬ 
lattice to qualify as a Weyl superconductor. Strong spin- 
orbit coupling is preferred because it gives a steep dis¬ 
persion around the node, making it well separated from 
other sub-gap excitations. 



FIG. 3. The dispersion of energy with respect to kx at the 
Weyl nodes for different values of spin-orbit coupling Vso- The 
other parameters are identical to the left panel of Fig. 2. 



FIG. 4. The pseudospin h (arrows) on equal energy contours 
(grey lines) near a pair of Weyl nodes. Note that ky = 0 and 
/lx = 0. To show the (three dimensional) h vector on the 
{kx, kz) plane, we plot hz along the kz axis and hy along the 
kx axis. Panel A is near the node at —k^ whereas panel B 
is near the node at A = 0.05/i, Vso = 0.5/i, Vz = 0.05/i, 
d/n = lO/kr- 



FIG. 5. The pseudospin h on equal energy contours near two 
pairs of Weyl nodes, G, D, E, and F. ky = 0 and hx = 0, hz 
is along the kz axis and hy is along the kx axis. A = 0.05/i, 
Vso = 0.5/i, Vz = 0.05/i, d/jv = Q.b/kp- 


given node. 



The three components of the h vector are the expectation 
values of the corresponding Pauli matrices fi in state |x), 

= ^{\u\^-\v^\), (13) 

h^=Re(?x*u), (14) 

h^=Im(i4*u). (15) 

The result is illustrated in Fig. 4 for the case of single 
pair of Weyl nodes. The pseudospin texture suggests that 

= V\\ [qxfy + qyfx] ± VzqzTz (16) 

for the two nodes at respectively. Using the formula 
for JF above, we can rescale q to find their correspond¬ 
ing topological charge ^ = ±1, i.e., the two Weyl nodes 
have opposite topological charge (chirality). The char¬ 
acteristic spin texture around two pairs of Weyl nodes 
is illustrated in Fig. 5. We observe that the topological 
charges of the two nodes on the positive kz axis, such 
as E and F, are opposite of each other. The two nodes 
at mirroring position e.g., D and F, also possess 

opposite charge. 


III. FERMI ARC 


We have numerically computed the pseudospin texture 
near the Weyl nodes. Consider a state |Tk) with positive 
energy and a momentum k on one of the cones shown in 
Fig. 2, and construct a spinor y by projecting |Tk) onto 
the |T±) basis formed by the zero energy states at the 


The Andreev bound states formed at the surface of a 
Weyl superconductor are very peculiar. The zero energy 
surface states in the two dimensional momentum space 
(such as the surface Brillouin zone) take the shape of a 
continuum line connecting two Weyl nodes of opposite 
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topological charge. Its topological origin has been dis¬ 
cussed in the context of Weyl semimetal^^’^^ and the A 
phase of 3p[022,23 repeated here. Instead 

we explicitly compute the spectrum of S-M superlattice 
structures that have a finite width L in the ^-direction 
to demonstrate the Fermi lines. Since the wave function 
has to vanish at ^ = 0 and y = L, we can expand it in 
sine Fourier series, 


y, z) = ^ sm{—y). (17) 

G,n 

Then the BdG equation becomes 

'Ho{G,n)^G,n + E '^nn' ^G'n' = e4>G,n (18) 

n',G' 


where we have defined 


Ho(G,n) 


^{kx,n,kz^G) Aia2 

-Aia2 -^{kx,n, kz + G) 


(19) 


with kz +^G) = {kl + + {kz + G)‘^)l2me - /i, 

kfi — TLTT j T, and Vnn' (p) — ^nn' [ ^so^xGy T '^zGz] G T 3 T 
Vsol^nn'O-x with 


e/A 
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2 

Knn' = I J [-idy] sin (k'^y). ( 20 ) 

In the numerical solution of the BdG equation, Eq. 18, 
it is important to keep enough kn- 

Fig. 6 shows the calculated finite-slab spectrum for 
single pair (upper panel) and two pairs (lower panel) of 
Weyl nodes. In the former case, a continuous line of 
zero energy states is formed along the kz axis connecting 
the bulk Weyl node at k^ to that at —k^. These Fermi 
arc (in this case just a straight line) states are absent in 
the bulk spectrum and correspond to the surface states 
at ^ = 0 and L. In fact, a very small gap is visible 
due to the hybridization of the two surfaces that are of 
finite distance L apart. In the latter case, the Fermi arc 
connects the two bulk Weyl nodes on the positive kz axis. 
(There is another arc along the negative kz axis, which 
is not shown.) This is in accordance with the topological 
charge of the nodes found above from the spin texture. 

IV. PHASE DIAGRAM 

The spin texture and the Fermi arc surface states pre¬ 
sented above unambiguously established the existence of 
pairs of Weyl nodes in S-M periodic structures. We have 
systematically scanned the parameter space of this model 
and the resultant phase diagram is shown in Fig. 7. Most 
strikingly, one observes a series of lobes (light gray re¬ 
gions, Wi) that feature single pair of Weyl nodes, sim¬ 
ilar to that shown in the left panel of Fig. 2. Roughly 
speaking, each lobe appears when kpd = nir {n E Z) 
and Vz > A. The dark black regions (W 2 ) represent 


FIG. 6. Fermi arcs on the surface of the S-M superlattice. 
The plot shows the low energy spectrum e{kx = 0,kz) of a 
finite slab y E [0,L]. Top panel: the arc connects kP and 
—kP (only the positive kz is shown), kpd — 50, Vz — 0.6A. 
Bottom panel: the Fermi arc along the kz axis connects two 
Weyl nodes on the positive kz axis, kpd — 30, Vz — 1.1 A. In 
both cases, kpL = 200, y = 20A, Vso = A. 


phases with two pairs of Weyl nodes. They also appear 
as a regular array but are well separated from each other 
and much smaller in area compared to Wi. In fact, W 2 
can be viewed as the overlapping regions of two adjacent 
Wi phases, as seen in the inset of Fig. 7 which shows 
the details near d/ir = 10. In the rest of the phase dia¬ 
gram (white regions. Wo), the spectrum is gapped, even 
through the gap may be numerically small (see for exam¬ 
ple the spectrum at point A shown in the first sub-panel). 

The evolution of the spectrum along a vertical cut in 
the phase diagram, from point A to point F, is illustrated 
in the sub-panels of Fig. 7. As Vz is increased at this par¬ 
ticular value of d, the lowest branch of the spectrum is 
pushed down by the increasing Zeeman splitting to touch 
E = 0, entering the W 2 phase (A^B). With further in¬ 
crease of Vz, the newly born twin of Weyl nodes with 
opposite charge become increasingly detached from each 
other, with one heading towards kz = 0 and the other to¬ 
wards the BZ boundary. The latter gets gapped out once 
it reaches the BZ boundary where it annihilates with its 
mirror image living on the negative kz axis, thus mark¬ 
ing the transition from the W 2 phase to the Wi phase 
(G^D). Further increase of Vz will push the spectrum 
away from E = 0, and the two Weyl nodes of opposite 
charge annihilate with each other at kz = 0 , leaving be- 
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the material-specific properties, especially regarding the 
coupling between S and M, and in the limit where the 
thickness of S and M is comparable to each other. As 
a result, the lattice model can potentially provide more 
quantitative predictions of the design parameters of Weyl 
superconductors. We will illustrate this approach by fo¬ 
cusing on the superconductor-magnetic topological insu¬ 
lator (S-TI) superlattice structure proposed in Ref. 17. In 
contrast to Ref. 17, however, we start from microscopic 
models of S and TI instead of the low energy surface 
degrees of freedom (i.e., the Dirac electrons). 

For simplicity, we model both S and TI on cubic lattice 
with lattice constant a. Each unit cell of the superlattice 
consists of Ns layers of S and Nt layers of TI stacked 
along the ^ direction. Let i be the layer index, the eigen¬ 
value problem has a tri-diagonal structure 

= {E- (21) 

Here is the hopping matrix coupling layer i to the 

neighboring layer i + 1, Hi is the Hamiltonian for the 
Tth layer, and is the wave function at the Tth layer. 
Note that the transverse momentum ky = {k^^ky) is con¬ 
served. For each S layer, i.e., i G [1, A^^], 

ff.(k„) = ffs(k,|) = ( ) , (22) 

where C(k||) = —2ts{coskx + cosky) — fig and k is mea¬ 
sured in units of 1/a. The hopping between two adjacent 
S layers is simply 


= ( ^o' t. ) ■ 

For example, we take tg = O.lSeV, and /i = —4tg. We 
consider Bi 2 Se 3 as a prime example of 3D Z 2 topological 
insulators, and model each TI layer by^^ 


FIG. 7. Top: Phase diagram of the nodal structure of the 
S-M superlattice as a function of superlattice spacing d and 
Zeeman field Vso in the regime of sufficiently strong spin-orbit 
coupling. The shaded (light grey) region corresponds to one 
pair of Weyl nodes, and the dark black region corresponds to 
two pairs of Weyl nodes. The spectrum is gapped in the white 
region. A = 0.05/i, the superlattice period d is measured in 
l/kp- Inset: zoom-in of the region around dkp — 10. The 
energy spectra (the lowest energy) along a vertical cut, from 
point A to point F, are shown in the bottom panel. 


^M(k||) = mro+a 2 sin/ca,ri+a 2 sin/c^r 2 +'i;;s(J 3 (g)i, (24) 

where m(k||) = M — 2bi -h 262 (cos-h cosky — 2 ). We 
choose the basis (|p+ f), |p+ |), |p_ t)/P- i))^ where p± 
labels the hybridized Pz orbital with even (odd) parity^^. 
The Gamma matrices are defined asFo =T 3 ( 8 ) 1 , F^ = 
with Ti {cFi) being the Pauli matrices in the orbital 
(spin) space. Vz is the Zeeman splitting for magnetically 
doped Bi 2 Se 3 ^^“^^. The coupling between two adjacent 
TI layers is given by 


hind a vacuum (E^F). 


tM — ^iTo — 


(25) 


V. LATTICE MODEL 

Now we turn to S-M heterostructures described by 
tight binding Hamiltonians defined on discrete lattices. 
Compared to the continuum model above, the lattice 
model can sometimes offer more realistic descriptions of 


The isotropic version of Km and with ai = a 2 , = 

62 , was proposed by Qi et al as a minimal model for 3D 
topological insulators^^. To mimic Bi 2 Se 3 , we set the 
lattice spacing a = 5.2A, which gives the correct unit 
cell volume, and a^ = Aija^ bi = Bija^ for i = I, 2. The 
numerical values of M, A^, Bi are given in Ref. 25. With 
these parameters, our model yields the correct band gap 
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and surface dispersion, it also reduces to the continuum 
k • p Hamiltonian (the Bernevig-Hughes-Zhang model) 
in the small k limit^^, aside from a topologically trivial 
eo(k) term. To describe the superconducting proximity 
effect, we have to generalize the TI hamiltonian above 
into the particle-hole space. For i G [Ns + 1, Ns + A^t], 



and accordingly. 




tM 

0 



(26) 


(27) 


Finally, the hopping from S to TI is a 4 x 8 matrix. 


' Ns,Ns+i 


_ ( isMO^w) 0 


V 0 


^Sm{ ^||) 


(28) 


with 


+ X _ / J+ 0 J_ 0 

tsM(k||)-l 0 J, 0 J_ 


(29) 


Here J± is the overlap integral between the p-orbital p± 
of TI and the s-\ike orbital of S. For simplicity, we assume 
the spin is conserved during the hopping, and from the 
orbital symmetry, J_^ = —J_ = J where J can be tuned 
from weak to strong^^. Small J mimics a large tunneling 
barrier between S and TI, while large J describes good 
contact, i.e., strong coupling between S and TI. 

A standard lattice Fourier transform from the layer 
index i to quasi momentum kz inside the reduced Bril- 
louin zone, [—7r/A^a, 7r/A^a] with N = Ns + gives 
each a phase factor The Hamiltonian of 

one unit cell of the superlattice is a matrix of the size 
ANs + 8A^t subject to periodic boundary conditions. 
Then a numerical diagonalization yields the band struc¬ 
ture E(kx^ ky^kz) of the S-TI superlattice. 


VI. ANDREEV BOUND STATES AT THE S-TI 
INTERFACE 

Before discussing the S-TI multilayer system, it is 
worthwhile to first consider the spectral properties of 
a single S-TI interface which have been studied exten¬ 
sively since the pioneer work of Fu and Kane^. Com¬ 
parison with these known results will serve as a critical 
check of our lattice model presented above. It also es¬ 
tablishes the connection between the microscopic model 
here and the effective model of Fu and Kane for the S-TI 
interface^, as well as that of Meng and Balents for the 
S-TI superlattice^^. 

The spectrum of a single S-TI interface can be conve¬ 
niently extracted from our lattice model by taking var¬ 
ious limits. Firstly, by setting J = 0 and Vz = ^ (but 
keeping finite A), the low energy spectrum reduces to the 


surface states of the topological insulator. As the thick¬ 
ness Nt is reduced, the linear Dirac spectrum acquires a 
gap due to the hybridization between two TI surfaces^^. 
For example, the gap is around 0.012eV when Nt = 8. 

E(eV) 



FIG. 8. Andreev bound states at the S-TI interface. Ns = 
Nt — ky — kz — The blue, green, and red curves are 
for J/ts — 0.2, 0.5, 1 respectively. Here = 0, A = 14meV, 
ts — O.lSeV, and p = —4ts. 


Secondly, by setting Vz = 0 and keeping Nt and Ns 
large, the low energy spectrum E{kx, ky^kz = 0) reduces 
to that of a single S-TI interface, since all the interfaces 
are sufficiently far apart and essentially decoupled from 
each other. Note that in this limit, the BZ is very small, 
and the dispersion along kz is negligible so we set kz = 
Fig. 8 shows the evolution of the subgap spectrum as J is 
increased from the tunneling to the strong coupling limit. 
In each case, the dispersion of the sub-gap modes can be 
fit well by formula ^(ky) = y^A^ + ('^s|k|| | E psY which 
follows from Fu and Kane’s phenomenological model^. 




idyNs \ 


(30) 


where hg describes the helical Dirac electrons^’^^. 


^s(k||) — fJjg T (81) 


This suggests that the Fu-Kane model is valid in a 
broad range of coupling strength between S and TL 
Yet, as clearly seen in Fig. 8, the effective parameters 
{ps, NsPs) ill Hfx depend sensitively on J. They may 
get strongly renormalized from their respective nominal 
values estimated from the bulk parameters by the prox¬ 
imity effect (the coupling to S). For device applications, 
e.g. for the generation and manipulation of Majorana 
zero modes, a large A^ and thus a strong S-TI coupling 
is preferred. In this limit, it is more natural to think of 
the interface state as the Andreev bound state which pen¬ 
etrates into the superconductor over the coherence length 
but decays rapidly (over a distance on the atomic scale) 
in the TI^^. The lattice calculation presented here is in 
agreement with the wave function calculation of Lababidi 
and Zhao in Ref. 20 and the Green function calculation 
of Grein et al in Ref. 31 . 









Thirdly, can be easily incorporated into the Fu-Kane 
Hamiltonian. For a single S-TI interface, it opens up a 
Zeeman gap at ky = 0. Using this as the starting point, 
Meng and Balents^^ analyzed the effective Hamiltonian 
of the S-TI superlattice and arrived at a very clean phase 
diagram in the plane of and A. 


VII. WEYL FERMIONS 

In search of Weyl nodes within our microscopic lat¬ 
tice model, we shall focus on “ideal” conditions provided 
that they seem experimentally feasible. High temper¬ 
ature superconductor (BSCCO) in proximity to Bi 2 Se 3 
was reported to induce a gap of 15meV, and the pair¬ 
ing symmetry was postulated to be s-wave because no 
d-wave nodes were observed^^’^^. In comparison, the gap 
of Bi 2 Se 3 grown on NbSe 3 is on the order of meV^^. We 
will consider a fairly large gap A = 14meV. The Zee- 
man field Vz can exceed the value of A. For example, 
chromium-doped Bi 2 (Sea,Tei_a ,)3 has an exchange gap 
of 40meV^^, and Bi 2 Se 3 doped with develops 

an exchange gap ranging from 10 to 60 meV. As to the 
number of layers for each materials, we will consider for 
example Ns = 5, Nt = 3 which gives a large BZ, al¬ 
lowing kz to have significant dispersion. Certain natural 
multilayer heterostructures (not superconducting) of TI 
such as (PbSe) 5 (Bi 2 Se 3 ) 3 ^ have been reported^^. Epi¬ 
taxial growth of TI films with 1 to 12 quintuple lay¬ 
ers on superconducting substrate has been successfully 
demonstrated^^. Also the superconductor BaBi 03 with 
Tc ^ 30K was predicted to turn into a TI upon electron 
doping^^. This implies that the S-TI superlattice may 
even be realized based on a single compound by modu¬ 
lated electric or chemical doping. 


E(eV) 



FIG. 9. One pair versus two pairs of Weyl nodes in the S- 
TI superlattice. The blue (orange) spectrum corresponds to 
Vz — 3A (5A). Ns — 5, Nt — 3, with J = 0.7ts, is — O.lSeV, 
A = 14meV. 

Fig. 9 shows the low energy part of E{kx = ky = O^kz) 
for Vz = 3A (blue) and 5A (orange). They have one pair 
and two pairs of Weyl nodes on the kz axis respectively. 
One can explicitly check that these nodes are the only 
zero energy states within the reduced BZ, and the energy 





FIG. 10. The number and location of Weyl nodes in the 
S-TI superlattice as functions of the Zeeman field Vz. The 
circle, diamond, triangle, and square are for Ns = 5,7,9,11 
respectively. Nt = 3 with J = 0.7ts, ts = O.lSeV, A = 
14meV. 


E is indeed linear in k — k^. Fig. 10 summarizes the lo¬ 
cation (and number) of the Weyl nodes on the positive kz 
axis as the Zeeman field Vz is increased. The phases and 
phase boundary can be easily read off from Fig. 10. Take 
Ns = 5 (the empty circle) for example, for Vz < 1.6A, 
the spectrum is gaped and the systems is in phase Wq. 
For Vz G [1.6A,4A], there is only one node on the posi¬ 
tive kz axis. And the node moves away from = 0 as 
Vz is increased. This is phase Wi. For Vz > 4A, a second 
nodal point appears at the BZ boundary kz = 7r/{Na). 
The system enters the W 2 phase. Overall, the evolution 
of the nodal structure here is similar to that of the con¬ 
tinuum model discussed above. Fig. 10 also compares 
the phases for increasing number of superconducting lay¬ 
ers. As a general trend, the critical Vz required for the 
Wi ^2 to appear increases with Ns. In the limit of large 
Ns, the dispersion along kz becomes very flat. 


VIII. CONCLUDING REMARKS 

We have presented two complementary approaches to 
model and compute the properties of S-M superlattice 
structures. The first approach is based on a simple con¬ 
tinuum model where M is described by periodic spin ac¬ 
tive potentials that are spatially thin compared to the 
thickness of the superconductor, and accordingly, the size 
of the unit cell can be larger than the superconducting 
coherence length. The second approach is based on a 
tight binding lattice model describing alternative layers 
of S and TI, both of which can be only of several layers 
thick with a tunable coupling strength between the two 
materials. In both models, we find phases that have one 
pair or two pairs of Weyl nodes. Together with previous 
results based on tunneling Hamiltonians^^, our study un¬ 
ambiguously establishes that (a) S-M periodic structures 
can behave as Weyl superconductors at low energies; and 
(b) for this to occur, neither topological insulators nor 
Dirac electrons are necessary, only the right combination 
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of spin-orbit coupling and Zeeman splitting are required. 
These observations generalize the proposals of realizing 
gapped topological superconductors featuring Majorana 
zero modes in one and two dimensions^^^ to gapless topo¬ 
logical superconductors in three dimensions using peri¬ 
odic structures of S and M. We hope that the theoretical 
analysis presented here can stimulate experimental work 
to explore these ideas. 

The emergence of Weyl fermions at low energies out of 
the vacuum of a conventional 5-wave superconductor is 
quite striking. While these low energy quasiparticles can 
be viewed as the Andreev bound states formed at the S- 
M interface dispersing with kz and crossing zero energy 
at isolated k points, they are located in k space near 
k = 0 (for kpd ^ 1 ), instead of being around the origi¬ 
nal Fermi surface |k| = kp^ In other words, the spectral 
weigh is transferred not only from high to low energy. 


but also from high to low momentum, by the presence 
of the periodic spin active potential. In this paper we 
have only considered one dimensional S-M superlattice 
structures. It is very likely that equally interesting phe¬ 
nomena may rise for “metamaterials” of superconductors 
and spin active materials that have more complicated pe¬ 
riodic structures in even higher dimensions. 
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